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Digital  Mammographic  Image  Compression:  Structure 
Lossless  Mammogram  Compression  and  Preliminary 
Lossy  Compression  By  Wavelet  Transformation 


Jun  Wang 
Dr.  H.K.  Huang 

1.  Introduction 


1.1  Digital  Mammography 

Breast  cancer  is  the  second  leading  cause  of  death  from  cancer  among  women  in  the  United 
States  1.  One  in  every  nine  women  will  have  breast  cancer  at  some  point  in  her  life  time.  Studies 
show  that  early  breast  cancer  detection  can  reduce  the  mortality  2.  Routine  screen-film  X-ray 
screening  has  been  recognized  as  the  best  means  to  detect  early  cancer  at  the  present  3.4,5 
Although  conventional  screen-film  mammography  has  been  shown  to  be  very  effective,  digital 
mammography  has  emerged  as  a  new  technology  that  will  overcome  some  of  the  limitations  of  the 
screen-film  mammography  and  offers  other  additional  benefits 

Digital  mammography  offers  the  following  advantages.  First,  screen-film  mammographic 
systems  produce  poor  contrast  images  for  dense  breasts  because  the  glandular  tissues  and  lesions 
are  essentially  the  same  density.  Digital  detectors  have  wider  dynamic  range  and  images  can  be 
post-processed  to  offer  even  more  contrast.  This  makes  direct  digital  mammography  an  attractive 
way  to  image  dense  breasts.  Second,  digital  image  processing  techniques  allow  the  change  of  the 
image  contrast,  enhance  lesion  detection,  and  optimize  image  display.  Such  techniques  thus 
provide  better  image  quality  to  support  higher  diagnostic  accuracy.  Third,  computer-aided 
diagnosis  (CAD)  involves  a  computerized  analysis  of  radiographs  that  can  be  used  as  a  "second 
opinion"  by  the  radiologist.  For  mammography,  CAD  offers  many  practical  schemes  to  detect  the 
masses  and  clustered  microcalcifications  on  mammograms  Fourth,  telemammography  provides 
a  fast  way  to  access  experts  over  long  distance  for  consultation  10. 


There  are  two  ways  to  generate  a  digital  mammogram:  a)  secondary  digitization,  in  which 
conventional  films  are  digitized;  and  b)  acquisition  of  primary  digital  images.  Previous  studies 
showed  that  at  least  100  p.m  sampling  distance  is  required  in  order  to  achieve  an  adequate 
resolution  for  digital  mammography  ^  ^12  Digitization  creates  9  Mbytes  of  data  at  the  resolution  of 
100pm  for  a  7'X9'  film,  and  36  Mbytes  of  data  at  a  resolution  of  50  pm.  Storage  cost  and  the 
transmission  cost  of  such  large  images  will  be  a  problem  as  more  digital  mammograms  are  created. 
Image  compression  reduces  the  data  storage  while  preserving  the  useful  information.  It  provides  a 
solution  to  handling  such  large  digital  images  efficiently. 

1.2  Image  Compression 

Image  compression  is  a  technique  that  reduces  the  number  of  bytes  needed  to  store  an  image. 
As  a  result  it  also  reduces  the  transmission  time  of  the  image.  The  goal  of  image  compression  is  to 
use  a  minimum  number  of  bytes  to  store  an  image  while  preserving  the  original  image  quality. 
There  are  two  types  of  image  compression:  lossy  and  lossless.  Lossless  compression  can  achieve 
a  compression  ratio  of  around  2  or  3  for  medical  images,  and  will  allow  exact  reconstruction  of  the 
original  image.  Lossy  compression  can  achieve  much  higher  compression  ratios,  although  the 
original  image  can  not  be  exactly  reconstructed. 

Various  lossless  and  lossy  image  compression  techniques  have  been  developed  in  the  past 
years  13-17.  Most  of  them,  however,  are  for  engineering  applications,  and  do  not  consider  the 
stringent  requirements  of  image  size  and  quality  in  medicine.  Although  several  investigators  have 
studied  compression  techniques  for  medical  images  18-22^  few  such  techniques  concern 
mammograms.  Kuduvalli  investigated  performance  of  several  reversible  image  compression 
techniques  for  high  resolution  teleradiology  23.  The  compression  schemes  implemented  are: 
Huffman  coding.  Arithmetic  Coding,  Lempel-Ziv  coding,  two-dimensional  linear  predictive 
coding,  and  transform  coding.  These  coding  techniques  are  applied  to  a  few  maiiunograms  and 
chest  radiographs  digitized  up  to  4k  X  4k,  10  bits/pixel.  The  results  show  that  multichannel  linear 
prediction  gave  the  best  results  at  an  average  pixel  rate  of  2.83  bits/pixel  for  five  mammograms. 


Benoit-Cattin  applied  a  wavelet  transform  to  mammograms  24.  Then  JPEG  was  applied  to  code 
only  the  low  resolution  subimage.  They  found  that  the  hybrid  method  achieves  a  lower  bit  rate 
than  solely  the  JPEG  method  at  the  same  SNR.  Lucier  evaluated  a  Harr  wavelet  transform  coding 
on  mammograms  25.  Images  with  clustered  microcalcifications  were  compressed  at  two  different 
ratios,  and  the  image  quality  was  evaluated.  The  results  show  that  wavelet  compression  could  be 
used  to  achieve  a  high  compression  rate  of  mammograms  without  losing  small  details  such  as 
microcalcifications. 

Mammograms  have  unique  characteristics  that  previous  work  did  not  consider:  a)breast  images 
have  a  regular  round  shape,  the  background  area  is  relatively  large  in  a  rectangular  film;  b)  breasts 
contain  mostly  soft  tissue,  so  image  contrast  is  less  compared  with  other  types  of  images.  In  this 
work,  we  develop  a  mammogram  compression  method  that  takes  advantage  of  these 
characteristics. 

1.3  The  Goal  Of  The  Present  Work 

The  purposes  of  this  research  are: 

1)  Develop  a  structure  lossless  compression  technique  for  mammograms  by  incorporating  the 
breast  imaging  shape  and  contrast  characteristics.  The  algorithm  developed  combines  segmentation 
with  a  lossless  compression  scheme  to  achieve  higher  compression  ratio  than  those  of  other  types 
of  lossless  compression  methods. 

2)  Implement  a  wavelet  lossy  compression  method  for  mammograms. 

1.4  The  Approach 

Structure  lossless  compression  has  three  steps.  The  first  step  identifies  the  orientation  of  the 
image,  i.e.  which  side  of  the  image  that  the  chest  wall  is  located.  The  second  step  segments  the 
breast  image  from  its  background.  The  third  step  compresses  the  portion  of  the  image  that  contains 
the  breast  tissue  with  the  lossless  compression  technique. 


2.  Method 


2.1  Digitization 

We  selected  14  cases  from  the  UCSF  mammography  teaching  file.  A  total  of  50  films  were 
digitized  at  the  resolution  of  100  iim  and  200  |J.m.  A  resolution  of  200  lam  may  not  good  enough 
for  clinical  purposes.  It  is  so  chosen  simply  because  images  digitized  at  this  resolution  have  small 
sizes  and  are  convenient  for  algorithm  development. 

A  sheet  of  film  contains  an  image  of  a  breast,  a  label  indicating  screening  position,  and 
possibly  patient  information.  Film  size  is  normally  7X9  inches  or  10  X  12  inches.  The  size  of  a 
breast  image  in  a  film  may  vary  a  great  deal.  An  image  may  only  occupy  a  small  portion  of  the 
film.  We  only  digitized  a  rectangular  area  that  contains  the  breast  image  rather  than  whole  film. 
The  patient  information  and  position  information,  which  are  located  on  the  edge  of  the  film,  are 
often  not  be  digitized.  This  information  can  be  input  as  a  text  file  from  the  image  header. 

The  average  image  size  is  1.3  Mbytes  at  200  |xm,  and  6  Mbytes  at  100  |xm. 

2.2  Orientation  and  Segmentation 

The  first  step  of  the  compression  process  is  to  find  the  orientation  of  the  image,  i.e.  which  side 
of  the  image  is  the  chest  wall.  Mammograms  normally  have  a  Cranial-Caudal  (CC)  view  and  a 
Medial  Lateral  Oblique  (MLO)  view  shown  in  Figure  1.  The  background  is  normally  dark  and  the 
portion  with  the  breast  tissue  is  normally  gray.  The  chest  wall  can  be  either  on  the  left  side  or  the 
right  side  of  the  image,  depending  on  the  way  that  the  film  is  digitized. 

To  determine  the  orientation  of  the  image,  ten  reference  points  along  the  left  and  the  right  edges 
of  the  image  was  selected  as  shown  in  Figure  2.  The  x  positions  of  the  left  and  right  reference 
points  are  located  on  5%  and  95%  of  total  width  of  the  image,  respectively.  The  y  positions  of  the 
five  reference  points  are  located  on  10%  to  90%  of  total  length  at  a  20%  increment.  The  average 
gray  level  of  each  reference  point  is  calculated  by  averaging  the  surrounding  3X3  pixels.  If  a 
reference  point  is  located  inside  the  breast,  the  gray  level  is  higher  than  if  a  reference  point  is 


located  in  the  background.  Comparing  the  gray  level  of  the  reference  points  on  the  left  and  right 
side,  we  can  determine  the  orientation  of  the  breast  image  and  the  background  gray  level. 


Left  CC  Right  CC  Left  MLO  Right  MLO 

(a)  (b) 

Figure  1.  Two  different  views  of  mammograms  (a)  CC  view  (b)  MLO  view 


Figure  2.  The  positions  of  the  reference  points 

The  next  step  of  the  compression  method  is  to  find  the  boundary  of  the  breast  tissue  and 
segment  the  breast  from  the  rest  of  the  background.  We  first  analyze  the  gray  level  distribution  of 
a  mammogram.  Figure  3  shows  a  profile  of  a  horizontal  scan  line  of  an  right  CC  image.  The  left 
side  is  the  chest  wall  and  the  right  side  is  the  background.  The  gray  level  of  the  image  is  around 
2000-3500  and  drops  to  a  background  of  around  200-300.  Near  the  edge  of  the  breast,  the  gray 
scale  decreases  gradually. 
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Figure  3.  One  scanning  line  of  an  Right  CC  mammographic  image 

We  use  a  thresholding  method  to  find  the  boundary.  The  edge  of  the  breast  is  not  clearly 
defined,  because  of  background  noise.  The  first  step  is  to  find  an  iso-gray  level  line  near  the  edge 
of  the  breast.  We  setup  a  thresholding  level  about  n  times  the  background  gray  level.  The  ratio  n 
is  adaptively  determined  by  the  ratio  of  the  average  image  gray  level  and  the  average  background 
gray  level  for  each  image.  The  image  is  scanned  line  by  line  to  find  the  location  that  reaches  the 
thresholding  level. 

The  first  scanning  line  is  started  from  one  third  of  the  height  of  the  image  from  the  top  instead 
of  starting  at  the  top  of  the  image  in  order  to  avoid  the  unpredicted  behavior  at  the  beginning  of  the 
image.  After  finding  the  first  thresholding  location,  the  scanning  continues  to  the  top  and  to  the 
bottom  of  the  image.  Since  the  breast  normally  has  a  smooth  boundary,  only  neighboring  pixels  of 
the  previous  thresholding  position  are  scanned  in  each  line. 

Since  the  threshold  position  is  not  the  true  boundary ,  the  second  step  is  to  expand  the  threshold 
position  further  to  the  background  region  until  the  true  boundary  is  reach  to  include  all  the  breast 
tissue.  Figure  4  shows  the  iso-gray  level  line  on  a  CC  view  image.  The  bold  line  is  the  true 
boundary  of  the  breast  and  the  thin  line  is  the  threshold  position  of  the  image  based  on  the 
thresholding.  The  size  of  the  expending  zoom  5  is  determined  by  the  following  procedure.  Ten 
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equally  spaced  horizontal  scan  lines  are  selected.  Each  scan  line  is  smoothed  and  the  distance 
between  the  threshold  point  and  the  edge  of  the  image  is  determined.  The  average  distance  is 
determined  from  the  10  scan  lines.  The  expending  zoom  6  which  equals  to  the  average  distance  is 

added  to  each  thresholding  position. 


Figure  4.  The  thresholding  of  a  CC  view  mammogram 


It  is  worth  mentioning  here  that  mammograms  normally  have  labels  on  films.  It  is  most  likely 
located  in  the  background.  When  we  segment  breasts  from  the  background,  the  labels  are  also 
removed.  In  a  digital  image,  we  can  overlay  text  on  the  image.  As  long  as  the  image  header 
contains  the  information  about  the  view,  it  does  not  affect  the  identification  of  the  image. 

2.3  Prediction  Compression 

The  next  step  is  to  compress  the  image  within  the  boundary  using  prediction  lossless 
compression.  X-ray  images  are  formed  due  to  the  different  attenuation  coefficients  of  tissue. 
Breasts  are  mainly  composed  of  adipose  and  fibroglandular  tissues.  Their  attenuation  coefficients 
differ  slightly  on  the  range  of  the  diagnostic  x-ray  energies.  Therefore,  mammograms  have  less 
sharp  edges  and  small  contrast  in  comparison  with  other  types  of  medical  images  as  breasts  contain 
primarily  soft  tissue.  Due  to  lack  of  edges,  we  select  prediction  as  an  average  of  the  previous 
pixels.  The  predicted  pixel  can  be  written  as, 

/(x,y)=  ^/(x-/,y-i)/« 

i,;>0 


(1) 


where  f{x,y)  is  the  predicted  pixel,  fix,y)  is  the  original  pixel,  i,  j  are  integers,  and  n  is  the 
number  of  the  pixel  summed.  The  prediction  error  can  be  expressed  as, 

e(x,y)  =  f(x,y)-f{x,y)  (2) 

e{x,y)  is  the  error.  After  the  prediction,  the  error  image  has  less  entropy  and  therefore  more 
efficient  for  the  entropy  coding.  Huffman  coding  is  selected  for  the  purpose  of  the  entropy  coding. 
Section  3.1  shows  the  results. 

2.4  Lossy  Compression  Method  With  Wavelet  Transform 

We  are  also  investigating  a  lossy  compression  scheme  using  wavelet  transform.  The  wavelet 
theory  has  been  developed  for  the  past  few  years  and  it  has  been  successfully  applied  to  image 
compression  26-30  a  wavelet  compression  method  normally  has  three  steps:  wavelet 
transformation,  quantization,  and  entropy  coding. 

The  basic  idea  of  wavelet  transformation  is  to  represent  any  arbitrary  function  as  a 
superposition  of  a  wavelet  basis  A  wavelet  transform  decomposes  a  signal  into  a  series  of 
smooth  signals  and  their  associated  detailed  signals  at  different  resolution  levels.  At  each  level,  the 
smooth  signal  and  associated  detailed  signal  have  all  the  information  necessary  to  reconstruct  the 
smooth  signal  at  the  next  higher  resolution  level.  The  transformed  signal  has  both  spatial  and 
frequency  information  from  the  original  signal  and  it  provides  a  good  representation  for  coding. 

The  quantization  step  reduces  the  data  precision  by  coarsely  quantizing  the  transformed  data. 
In  the  last  step,  entropy  coding  such  as  Huffman  or  arithmetic  coding  can  be  used  to  reduce  the 
coding  redundancy. 

2.4.1.  Wavelet  Transform 

The  wavelet  transform  can  be  implemented  by  a  two-channel  perfect  reconstruction  (PR)  filter 
bank  31.  A  filter  bank  is  a  set  of  filters,  which  are  connected  by  sampling  operators.  Figure  5 
shows  an  example  of  a  two-channel  filter  bank  applied  to  a  one-dimensional  signal.  x{n)  is  an 
input  signal.  Hg  and  are  analysis  filters  and  Go  and  G,  are  synthesis  filters.  is  a  lowpass 


filter  and  //,  is  a  highpass  filter.  If  the  output  x{n)  =  x(n-l),  where  I  is  a  delay,  then  the  two- 
channel  filter  bank  is  called  a  perfect  reconstruction  (PR)  filter  bank. 


x(n) 


1  Ho  H'l'  2  1 - 

1-2  -[^ 

1  Hi  |— 2  1 — 

— FnHGr 

x(n) 


2  Downsampling 


'[  2 


Upsampling 


Figure  5.  The  perfect  reconstruction  filter 

The  wavelet  transform  of  a  signal  can  be  obtained  by  repeatedly  applying  a  PR  filter  bank  to  the 
signal  in  a  pyramidal  scheme.  In  a  two-dimensional  case,  a  two  level  wavelet  decomposition 
process  is  shown  in  Figure  6: 


X  direction 


y  direction  x  direction 

keeping  every  other  pixels 


y  direction 


Figure  6.  Two  dimensional  wavelet  decomposition 

We  first  convolute  each  line  in  the  x-direction  with  filters  and  respectively,  followed  by 
subsampling  every  other  pixel.  We  then  convolute  the  resulting  signals  with  Hq  and  H^  in  the  y- 
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direction,  followed  by  subsampling.  Because  Hq  is  a  lowpass  filter,  and  is  a  highpass  filter 
the  resulting  signals  are  a  smooth  signal  ,  and  detailed  signals  /  ,  respectively.  The 

smooth  signal  /„+,  is  again  convoluted  with  filters  and  H,  respectively  in  x  and  y  directions, 
sampled  at  every  other  pixel,  creating  and  f' ^^2^.  The  same  process  can  be  continued  until 
the  desired  level  is  reached.  The  result  of  the  two  level  decomposition  contains  and 

f'  s  f  o  is  the  smooth  signal  and  f  and  /„+iS  are  detailed  signals  in  the  different 

frequency  bands. 

Figure  7  shows  one  mammogram  and  its  one  level  and  two  levels  of  the  wavelet 
decomposition.  Figure  7(a)  is  the  original  image.  After  the  first  level  wavelet  decomposition 
(Figure  7(b))  the  image  contains  four  components.  The  upper  left  comer  is  a  smooth  image,  the 
rest  of  the  comer  images  are  detailed  images.  If  the  image  has  sharp  edges  on  the  x  or  y  direction, 
the  upper  right  comer  will  contain  the  edge  in  the  y  direction,  the  lower  left  comer  will  contain  the 
edge  in  the  x  direction,  and  the  lower  left  will  contain  the  edge  in  the  diagonal  direction.  Since  the 
mammogram  does  not  contain  sharp  edges.  The  detailed  images  do  not  show  the  visible  edges. 
The  upper  right  comer  image  can  be  further  decomposed  similarly  into  four  more  components  at 
the  next  level  of  resolution  (Figure  7(c)). 

In  general,  the  low  frequency  component  of  the  wavelet  transform  is  about  l/2^"  of  the 
original  image  size,  but  contains  more  than  99  %  of  the  total  energy,  where  M  is  the  level  of  the 
decomposition.  The  high  frequency  components  are  separated  into  different  resolution  levels.  At  a 
particular  resolution  level,  each  block  contains  the  high  frequency  information  m  certain  directions. 
Blocks  at  different  levels  contain  similar  structure  but  different  frequency  band  information. 

2.4.2.  Quantization  and  Entropy  Coding 

The  second  step  of  the  3D  wavelet  compression  is  quantization.  The  purpose  of  quantization  is 
to  reduce  data  entropy  by  compromising  the  precision  of  the  data.  Reducing  entropy  allows  more 
compression.  The  quantization  step  maps  a  large  number  of  input  values  into  a  smaller  set  of 
output  values.  The  original  data  caimot  be  recovered  exactly  after  quantization.  It  is  therefore  very 
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Figure  7.  A  maminograin  and  its  wavelet  transformation,  (a)  Original  image,  (b)  One  level 


important  to  design  a  quantization  strategy  which  selectively  quantizes  the  wavelet  coefficients  and 
preserves  the  image  quality.  A  simple  uniform  quantization  is  implemented  in  this  study. 

In  the  third  step,  run-length  coding  followed  by  Huffman  coding  is  applied  to  the  quantized 
data.  Run-length  coding  is  effective  when  there  is  more  than  one  pixel  with  the  same  gray  level  in  a 
sequence.  This  method  uses  two  integers  to  represent  a  sequence  of  the  same  gray  level.  The  first 
integer  represents  the  length  of  the  sequence,  and  the  second  integer  represents  the  gray  level  of  the 
sequence.  The  longer  is  a  sequence,  the  more  efficient  the  run-length  coding.  Since  thresholding 
of  the  high  resolution  components  results  in  a  large  number  of  zeroes,  run-length  can  be  expected 
to  significantly  reduce  the  data. 

3.  Result 

3.1  Structure  Lossless  Compression 

The  structure  lossless  compression  method  was  applied  to  all  100  digitized  mammograms  at  a 
resolution  of  100  |im  and  200  |im.  The  segmentation  algorithm  correctly  segmented  breast  images 
from  their  backgrounds  for  each  of  the  100  mammograms.  We  applied  five  different  predictions. 
The  predictions  are  one  of  the  previous  neighboring  pixels  a,  b  ,  c,  or  their  combinations 
{a  +  c)/2,  (a  +  b  +  c)/ las  shown  in  Figure  8. 


b  c 


X  X 

a 

X  0 

Figure  8.  The  distribution  of  the  prediction  pixels 

Table  1  summarizes  the  results.  The  compression  ratios  are  higher  at  100  |im  than  those  at  200 
|im  for  all  the  predictors.  Column  5  gives  the  best  compression  ratios.  Comparing  the  first  three 
columns,  predictor  a  gives  the  best  results.  This  may  be  the  partially  caused  by  digitization 


process.  This  is  because  the  digitizer  scans  an  image  line  by  line,  an  immediate  previous  pixel  of 
the  same  scanning  line  has  better  correlation  with  the  current  pixel. 


Predictions 

a 

b 

c 

(a  +  c)/2 

(a  +  b  +  c)/2 

200  |im 

3.17 

3.07 

3.01 

3.24 

3.19 

100  |im 

3.34 

3.24 

3.18 

3.42 

3.37 

Table  1.  Summary  the  average  compression  ratios  for  different  predictors 


3.2  Wavelet  Compression 

We  also  applied  the  wavelet  compression  to  a  mammogram  at  different  compression  ratios. 
The  energy  distribution  of  the  wavelet  transformed  image  was  calculated.  The  total  energy  of  a 
digital  image  is  defined  as  the  summation  of  the  scjuare  of  the  pixels.  Table  2  lists  the  energy 
distributions  at  the  different  levels  of  decomposition. 


Low  Resolution 

Hiqh  Resolution 

level  1 

level  2 

level  3 

level  4 

1  level 

99.9844 

0.0156 

2  level 

99.9750 

0.0155 

0.0095 

3  level 

99.9629 

0.0154 

0.0093 

0.0124 

4  level 

99.9343 

0.0151 

0.0092 

0.0122 

0.0292 

Table  2.  The  energy  distribution  of  the.  wavelet  transformed  image. 


The  first  column  shows  the  number  level  of  the  decomposition.  The  second  column  shows  the 
percentage  energy  in  the  low  resolution  component  in  that  level.  The  rest  columns  show  the 
percentage  energy  at  different  levels  of  the  decomposition.  The  energy  in  the  low  resolution  area  is 
more  than  99.9  percent  for  all  level  of  the  decomposition  and  the  total  energy  in  the  high  resolution 
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areas  is  less  than  0.06  percent.  Therefore,  selectively  quantizing  the  high  resolution  components 
will  only  introduce  small  errors. 

The  decompressed  image  quality  is  measured  by  the  Peak  Signal  to  Noise  Ratio  (PSNR) 
defined  as: 


PSNR  =  20log 


where  maximum  gray  level  of  the  image  set,  N  is  total  number  of  pixels,  f(x,y)  is  the 

original  image,  and  f^ix^y)  is  the  decompressed  image.  The  denominator  is  the  Root  Mean  Square 
Error  of  the  decompressed  image.  The  larger  the  PSNR,  the  better  the  decompressed  image  quality 
is. 

The  compression  performance  is  measured  by  Compression  Ratio  (CR)  which  is  defined  as 
Original  bits  /  pixel 

Compressed  bits  /  pixel  (4) 


Figure  9  shows  compression  ratio  versus  PSNR  at  different  wavelet  decomposition  levels. 
The  level  of  decomposition  ranges  from  1  to  4.  At  each  level  of  decomposition,  a  range  of 
compression  ratios  were  obtained  by  varying  the  quantization  numbers.  The  PSNR  is  used  to 
measure  the  decompressed  image  quality.  At  the  same  PSNR,  the  compression  ratio  increases  as 
the  level  of  the  wavelet  decomposition  increases.  This  is  because  at  higher  levels  the  image  energy 
is  concentrated  in  a  smaller  low  resolution  area,  and  the  high  resolution  components  are  separated 
into  different  frequency  bands,  it  is  therefore  convenient  to  quantize  the  high  resolution 
components.  The  compression  performance  does  not  increase  when  the  level  of  the  wavelet 
decomposition  is  higher  than  four. 
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Figure  9.  The  compression  ratio  versus  PSNR  at  different  wavelet  decomposition  level 

4.  Conclusions 

During  the  past  year,  we  investigated  both  lossless  and  lossy  compression  methods.  In  the 
lossless  compression,  we  developed  a  structure  lossless  compression  method  for  mammogram 
image  compression.  The  algorithm  utilizes  the  unique  shape  characteristics  and  image 
characteristics  of  mammograms.  The  combination  of  the  segmentation  and  the  prediction  coding 
enable  us  to  achieve  high  compression  ratios  and  without  losing  any  useful  information. 

In  the  lossy  compression,  we  developed  a  wavelet  method  which  formulates  the  background 
for  improving  image  quality  at  high  compression  ratios  for  mammograms.  Currently,  we 
implemented  a  limited  number  of  wavelet  filters  for  the  transformation  and  utilized  simple  uniform 
quantization.  In  the  next  stage,  we  will  study  the  wavelet  filter  characteristics  and  select  the  best 
filter  for  mammogram  compression.  Uniform  quantization  will  be  replaced  by  an  adaptive 
quantization  method  based  on  mammogram  characteristics. 

PSNR  and  compression  ratio  have  been  used  for  assessing  image  quality  and  compression 
performance.  However,  PSNR  cannot  represent  the  human  visual  system.  In  the  next  stage,  we 
will  propose  quantitative  parameters  for  better  assessing  image  quality. 
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